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Suppression of electron current AI through a ID channel of length L connecting two Fermi liquid 
reservoirs is studied taking into account the Umklapp interaction induced by a periodic potential. 
This interaction opens band gaps at the integer fillings and Hubbard gaps 2m at some rational fillings 
in the infinite wire: L —f oo. In the perturbative regime where m ^ Vc/L {vc : charge velocity), and 
for small deviations 5n of the electron density from its commensurate values —AI/V can diverge 
with some exponent as voltage or temperature V,T decreases above Ec = max{vc/L,VcSn), while 
it goes to zero below E^. This results in a non-monotonous behavior of the conductance. In the 
case when the Umklapp interaction creates a large Mott-Hubbard gap 2m S> Tl inside the wire, 
the transport is suppressed near half-filling everywhere inside the gap except for an exponentially 
small region ofV,T< TLexp{-2m/TL). 
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I. INTRODUCTION 



Recent developements in the nano-fabrication tech- 
nique have made the ID interacting electron systems an 
experimental reality, and its quantum transport proper- 
ties have been the subject of extensive studies both ex- 
perimentally [1-3] and theoretically [4-16]. In realistic ex- 
perimental set-ups, the quantum wire is attached to two- 
dimensional regions called reservoirs or leads. A metal- 
lic phase of the infinite wire is known as a Tomonaga- 
Luttinger liquid [|l^. To describe ID transport phenom- 
ena in the realistic configuration a model was recently 
formulated of the inhomogeneous Tomonaga-Luttinger 
liquid (ITTL) |,|,§. It recovers the conductance G = 
2e^/h observed experimentally even in the presence of 
the electron-electron interaction in the wire [|l| (below 
we use the units, where e = h = I), although the previ- 
ous calculations on an infinitely longwire [Q,^ predicted 
the renormalized conductance ( see ^ for the further de- 
velopment). Calculation of the conductance suppressed 
by a weak random impurity potential in this model [|ic| ] 
had agreed with both the previous theoretical prediction 
[H and an experiment [Q . 

In this paper we consider effect of openin g a spectral 
gap in the wire on the ID transport. Theory pT| ] predicts 
that in ID besides band gaps produced by a periodical 
potential in the wire at the integer fillings, a repulsive 
interaction between electrons opens Mott-Hubbard gaps 
at some rational fillings. Furthermore recently Tarucha 
et al. [|l^ succeeded in introducing the ID periodic po- 
tential with a periodicity of order 40nm into the wire 
2/im in length and 50nm wide. This induces Umklapp 
scattering. The electron density n can be continuously 
controlled by the gate voltage, and one can satisfy the 
half-filling condition within an accessible value of n. If 
this condition is satisfied the system will becomes a ID 
(doped) Mott insulator with the Mott-Hubbard gap m 
for the infinite length wire. Then it will offer an ide- 



alistic system to study the quantum transport in Mott 
and doped Mott insulators in ID. Earlier a similar ex- 
periment had been carried out by Kouwenhoven at al [^ . 
Their wire was 3/im in length, 250-400nm wide and the 
period of the structure 200nm. To reduce the impurity 
backscattering they put the set-up in a strong magnetic 
field and observed an interference structure at the inte- 
ger edge-state plateaus in conductance vs. gate voltage. 
This structure partly had been explained as a band gap 
appearance. 

Inspired by these works, we study theoretically the 

1 — V characteristic of the wire of length L connected 
to leads at different temperatures T taking into account 
the Umklapp electron-electron interactions, which can 
be directly compared with the experiments. We for- 
mulate our model in the section II and check it on the 
description of the band gap effect at different temper- 
atures. We calculate |jlq] the suppression of the cur- 
rent AI perturbatively in the Umklapp scatterings in the 
section HI. The results are summarized in Figs. 1 and 

2 for a threshold structure near half filling. There are 
two energy scales, i.e., the finite size energy Tl — Vc/L 
(vc : charge velocity) and Ethr ^ VcSn {Sn : the devia- 
tion of the electron density from its value at the filling 
v equal to 1/2) measuring the incommensurability. For 
T,V > Ec = max(vc/ L, Ethr), the suppression of the 
current —AI/V diverges as max{T, V)'^^-^ jf t}jg short 
range interaction constant g for the forward scattering is 
less than 3/4. For T,V < Ec, on the other hand, the 
suppression —AI/V goes to zero as T,V ^ Q. Then we 
predict the non- monotonous temperature and/or voltage 
dependence of /, which is the clear signature of the Umk- 
lapp scattering effect. For small values of g, expected 
when the screening length of the interaction deter- 
mined by the close metallic gates is much larger than the 
width of the channel d and g (x 1/ yJ\nXcJd II4 ] , we pre- 
dict a few more threshold singularities. These features 
could be observed experimentally by changing the gate 
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voltage, bias voltage, and temperature. 

In the next section IV, motivated by a recent discus- 
sion [|l9|,^ that in a non-perturbative case of the large 
Mott-Hubbard gap 2m ^ v^/L = Tj^ the conductance 
is strongly suppressed at any low energy (jm) inside the 
gap similar to the band gap case, we consider transport 
through a ID Mott-Hubbard insulator of a finite length 
L beyond perturbative approach. Our presentation fol- 
lows ||2^ where we used a special value of the low energy 
constant of the interaction to map the problem onto the 
exactly solvable models. We find current vs. voltage V 
at high temperature T > max{m, Tl) and at low energy 
T,V <Ti^. The result shows that for the strong interac- 
tion creating a large Mott-Hubbard gap 2m ^ Tl inside 
the wire, the transport is suppressed near half-filling ev- 
erywhere inside the gap except for an exponentially small 
region of T^,r < TLexp(-2m/TL) 

II. MODEL 

Our model can be derived following |^ from a 1 chan- 
nel electron Hamiltonian 

■H= fdc{J2ij+{x){-^~EF)M^) 
J ^-^ 2m* 

a 

^•^{x)p^{x) -\- \yimp{x) + Vj,eriod{x)\p{x)} (1) 

with the periodic potential V-periodix) ( period a) assumed 
to be weak enough to justify the perturbative consider- 
ation of the Umklapp backscatterings. The Fermi mo- 
mentum kp and the Fermi energy Ep is determined by 
the filling factor v s& v = kpa/ir and Ep ~ vpkp. 
In Eq. the function ip{x) = const x 9{x)9{L — x) 
switches on the electron-electron interaction inside the 



wire confined in < a; < L. This interaction is assumed 
to be local, as the close metallic gate used in experi- 
ments to form the wire inevitably screens the long range 
Coulomb. Contribution of the random impurity poten- 
tial Vimv{x)p{x) to the conductance has been considered 
in some results of which we will use below. Fol- 

lowing Haldane's generalized bosonization procedure p^ ] 
to account for the nonlinear dispersion one has to write 
the fermionic fields as ipaix) = yJkp/{2TT) ^ exp{i{n + 
l){kpx + (f)a{x)/2) + i9„{x)/2} and the electron density 
fluctuations as p{x) — J^P'yi^)^ Pa-i^) — [dxfpaix) + 
2kp]/{2'K)'^ exp{in{kpx + (l)a-{x)/2)} where summation 
runs over even n and 0(j,0cr are mutually conjugated 
bosonic fields [(j)a{x), Oa-{y)] = i2i:sgn{x — y). 

After substitution of these expressions into (|l|) and 
introduction of the charge and spin bosonic fields as 
4'c,s = [<P] i %/2 the Hamiltonian takes its bose-form 
Ti. — Tio + T^bs ■ Here the free electron movement modi- 
fied by the forward scattering interaction is described by 

with gdx) — g for x ^ [0, L] ( 5 is less than 1 for the repul- 
sive interaction and it will be assumed below ), gdx) = 1, 
otherwise and Vc{x) = vp/gdx). The constants in the 
spin channel gs = l,Vs = vp are fixed by SU{2) symme- 
try. Keeping only the most slowly decaying terms among 
others with the same transferred momentum one could 
write the backscattering interaction as 



E%, r m(j)dx) ^ v-^ i-T /4'six)-. UKpdx) ' 
/ dx\ > Um cos(2 kmprnx ~\ ^ — )+ > t/m cost — z^]cos(2kmpnix -\ r= — 

-'O V2 ,,r^^n V2 V2 ' 

" even m>0 ' odd m>0 



A difference in the transferred momentums 2kmF is 
brought by the periodical potential with the period a: 
kmF = kp — Til/ {ma)) ^ where I is an integer chosen to 
minimize | kmP \ ■ We have omitted the m = term of the 
first sum: It can contain only the spin field and cannot 
affect the current in the lowest perturbative order. The 



most singular is the m ~ 1 term of the second sum re- 
sponsible for opening the band gap in the infinite wire at 
i> ~ integer. The dimensionless coefficients [/,„ originate 
from (p{x). 

In the spinless case we should put (f>c = (j>a in @ and 
change the backscattering interaction to 



3l 

Vp 



I dx^U, 



1 C0s{2kmF'niX + UKpa-ix)) 



(4) 



To generalize our perturbative results of the part HI to 
the spinless case one just needs to transform m, kmF to 
■\/2m, kmF 1^/2 in the expressions written in the spin case 
for even m and take m arbitrary integer. 

It is instructive first to examine how the band gap 



shows out in transport properties of the wire filled with 
the non-interacting electrons. The model is equivalent to 
a Dirac equation with the mass switched on inside the 
wire: 
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(5) 



a=R,L=± 



with m — TTEpUi/2 in the spin case of Eq. ^ and 
m — irEpUi in the spinless case of Eq. (^. Here ipuiL) 
describes the right (left) chiral electron and fc^ is counted 
from nl/a at v I. The transport is determined by the 
transmittance D which is a function of the electron en- 
ergy £ counted from the one of the middle of the gap 



D{e) 



(6) 



where the analytical extension is assumed at e < m and 
t]^ = L/vp. In particular, the linear bias conductance G 
is equal to 



G(r,/io) = 



1 



2ttT 



de 



D{e) 



1 + cosh((£ - Mo)/T) 



(7) 



for spin electrons where /io = kpvp gives a deviation of 
the chemical potential of the wire from that of the ly — I 
filling. The conductance has two regimes of behavior. 
1 . High temperatures T > T^^m - The asymptotics to 
can be written as 



TT V 2Tl + cosh(^o/r) 



(8) 



after noticing that is a quickly oscillating function 
against the slowly varying temperature factor. The co- 
efficient B changes slowly from to ~ tt and is defined 
below in Eq. (p^). The expression (||) shows that at the 
high temperature the band gap produces a smooth well 



in the conductance vs. chemical potential which becomes 
more narrow and deeper as T decreases. 

2. Low temperatures T < T^- The conductance is 
about D{pL^)/'K. It is suppressed in the middle of the gap 

G(0) = [7r(l + (mt^)^)] and approach its maximum 
I/tt away from the gap oscillating with period Tr/i^. In 
the limit mtL 1 this interference structure is perfectly 
periodical, as 



AGi = -{mtL)' 



, sin^((fci? — 7rl/a)L) 
{{kp - T:l/a)LY 



(9) 



Therefore the conductance has L/a humps in between 
two neighbor band gaps. We expect this being correct at 
any finite m from comparison with a tight binding model. 
In our model it holds on asymptotically if T:Vc/{am) ^ 1. 
Then the model works well. 



III. NARROW GAP: PERTURBATIVE 
APPROACH 

In this section we assume that the dimensionless co- 
efficients Ura in Eq.(^) are small enough to justify per- 
turbative calculation of the current. The variation of the 
current due to the backscattering is given by : ZV = 
-t/(2V27r) / (k[dJc{x),n] =V2j (k{S/Sq>,{x))nbs- At 
finite voltage V applied symmetrically to neglect the mo- 
mentum transfer variation, the average of AI decomposes 
into sum of the different backscattering mechanism con- 
tributions < Aim > in the lowest order. The even m 
terms involving only (j)c field are equal to 



< Mm > = 



m 



' VP ' 



dt 



dxidx2 < e 



(10) 



The current operator Aim has a high energy scaling di- 
mension m?g/2 and a free electron {g = 1) behavior at 
low energy. We will see below that the integral (|lo| ) 
scales at low energy with (vn? — 1) exponent and with 
{vn?g — 2) exponent at high energy. The most singular 
behavior is due to Umklapp backscattering at m = 2 
with the threshold voltage V — Ethr — 2k2FVc going to 
zero at the half filling. We assume V, Ethr > below. 
Less singular correction with m = 4 could become rel- 
evant at the one and three quarters fillings and so on. 
Expressions for the odd m terms include additionally a 
spin field correlator < e^'l>s{xut)/V2^~t^,ix2,a)/V2 > ^^^^^ 
the integrals in (Il0|). The high energy dimension of Aim 
in this case is m^/2 + 1/2. The most singular behavior 



occurs to the m = 3 term at the one and two thirds fill- 
ings. It has two threshold energies Ethrc.s = 2k'ipVc,s for 
Vc Vg. Neglecting a change of the TLL compressibility 
produced by the Umklapp scattering we can relate |23] 
the threshold energy to a deviation of the chemical po- 
tential of the wire fj, at V = from that of the rational 
filling. Since Apc — 2kp/'!T and d^/dpc — TTVc/{2g) we 
gather Ethr.c = 2gfi. However, in an experiment it is the 
average of the electrochemical potentials of the leads but 
not /i that is known. The latter is proportional to ^ with 
the coefficient [1 -I- 2ge'^ / {irVcTiCg)] if the gate voltage is 
fixed 1^^. This coefficient is about 1 if the density Cg of 
the capacitance between the wire and the screening gate 
is large. 

Correlator of the charge field exponents e^'t>':i^-*-) ^ 
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evolution of which is specified by Tio, could be (air / j3Y^ / (Y\ . sinh^(7r(a:/f ± (i — ia))//3))) in the fol- 
compiled from the correlators of the uniform TL lowing way ||4[|l0|] 
liquid K{x,t) = K{x,t,g,Vc) {K{x,t,g,v) = 



< gi0.(^,t)g-^0c(^/,O) >^ ^(2; - 2/, i) ]J 



K{2nL,0) 



^^i^'Ki2nL±\x-y\,t)^ 
^ ,K{2{nL + a;), Q)K{2{nL + y), 0) . -r^"+i /2 ^ ,K{2{nL ~ x), Q)K{2{nL - y), 0) ^-r^— 1/2 



(11) 



n=0 



K^{2nL-x-y,t) 



Here /3 is inverse temperature l/T and a = 1/Ep is 
the ultraviolet cut-off. This complicated form comes 
about through a multiple scattering at the points of joint 
X = 0, L. As a result of the scattering the correlator 
< 4>c{x,t)(l)c(i),Q) > becomes an infinite sum of the uni- 
form correlators taken along the different paths connect- 
ing points X and y and undergoing reflections from the 
boundaries at a; = 0, L. Each reflection brings addi- 
tional factor r = (1 — g)/{l + g)- The similar correlator 



< e 



> for spin field is K{x — y,t^ Ill's)- 
Below we analyze the current corrections (|lo|) for high 
(T > l/tL ^ Tl) and low (T < T^) temperatures, re- 



spectively. 

l.High temperatures T > - The uniform correlator 
K{x,t) goes down exponentially if distance between the 
points I X I exceeds the inverse temperature. There- 
fore only paths with length less than /3 contribute to 
the correlator ([Tl|). This means that the high temper- 
ature form of the correlator ( pi] ) reduces to the first mul- 
tipher K{x — y,t) up to a factor (1 + 0{exp{—tL/ (5)). 
Neglecting 0(Tl/T) quantity we can extend integration 
over xi — X2 in ( |lO|) from —00 to -l-oo. Then calcula- 
tion of the m — 2 contribution reduces to finding Fourier 
transformation F2g{q, e) of the correlator K^{x, t, g, Vc)'- 



1 fUoE' 



2^F 

9 



tLj2TF2g{2Eu,r,±2V) ^ ^2{ 



± 



r{2g)g > Tl'-Ep' 



±,± 



2ttT 



(12) 



One can easily see its behavior making use of the following asymtpotics 
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< ZV2 >cx 



< A/, > 



2 >« 



2^2 



9 J n \ {{V + Etnr)T/Elfs-\ 



V » Ethr: T 
V « Ethr » T 

Ethr » V, T 
V , Ethr « T 



(13) 




1<T E,h, ~2E„, V 

FIG. 1. Scliematic voltage dependence of tlie high temper- 
ature current corrections produced by the m = 2 Umklapp in- 
teraction A/2. Solid lines, Ethr = 0; dashed lines, Ethr 2> T. 



ity in the current voltage dependence diverges as {{V — 
Ethr)/Tf^~^ if g < 1/2 (Figjl). It becomes stronger 
in the differential conductance dependence. At g > 1/2 
the differential conductance correction dG'2 behaves as 
-{V/Epf^-^, and saturates at -{T/Epf^^^ below T, 
if Ethr < T] otherwise, the correction shows diver- 
gence — {[V — Ethr)/T)^^~^ smeared over T scale near 
the threshold and becomes suppressed exponentially as 
—exp{—Ethr/T) below it. Eq. (|l^) predicts also that 
the linear bias conductance G as a function of Ethr pro- 
portional to fi and the gate voltage has a T wide well 
at half filling similar to the band gap case discussed in 
the section IL However, the depth of the well oc 
increases with decrease of T only if the repulsive inter- 
action is strong enough. This is a bit unexpected result 
since we know [ pT| that any repulsive interaction opens 
the spectral gap at half-filling in the infinite wire. It will 
be discussed further in the next section. 



These asymptotics show that the threshold singular- 



Generalization to the other even m expressions for the 
backscattering current of the spin electrons needs just 
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changing: Ag ^ m'^g ,Ak2F^2V 'm{2k„iF±V). In the 
case of the spinless electrons: Ag — > 2m^g , Ak2F ± 2V^ — > 
m{2kmF i V) with arbitrary integer m. In particular, we 
find again ^-dependence for decrease of the conduc- 
tance produced by the band gap opening in the spectrum 
of free electrons. The edge singularity is characterized by 



a half of the scaling dimension for A/,„ since only one 
chiral component of the field (j)c contributes. 

As to the odd m terms, the two threshold energies 
Ethrc.s = '^kmpVc.s bccomc distinguishable if their dif- 
ference exceeds T. The leading high-temperature current 
correction reads as: 



\{mEthrs - qvs,±mV - s)F^{qVc,e) (14) 



Substituting zero temperature form of Fa func- 
tion Fa{q,w) = 8[sin{na)r{l - a)]2(a/2)2° J]^(w ± 
q)'^^^9{wzLq) in jl^ ) one can gather that the current cor- 
rection behaves as —{V/Ef)™' at large voltage V > 
Ethrcs, has a leading singularity -{{V -Ethrc)/T)'^''s/'^ 
smeared over T scale near the first threshold and —{{V — 
Ethr s)/T)™- near the second threshold (we assume 

Vc > Vs). Below the lowest threshold it becomes expo- 
nentially suppressed. These singularities result in the di- 
vergences of the differential conductance or higher deriva- 
tives of the current in voltage. The threshold behavior of 
the TO = 3 term of the differential conductance correction 
is divergent at Ethr c if < 2/9 and at Ethr s if 5 < 1/6. 
On the other hand, it is unlikely that the splitting of the 



threshold energies could be observed at the integer fill- 
ing factor where opening the band gap eventually makes 
electrons non-interacting. 

2. Low temperatures T,V << Tl - With lowering tem- 
perature we should expect that above current correction 
dependencies will be modulated by a ttTl quasiperiodical 
interference structure ||l^,0 and also a new low energy 
scaling behavior of the current correction operators will 
appear at V, T < T^. The dominant contribution to the 
integral of (|lO| ) comes from long times t^ti,. One can 
neglect the spacious dependence compared with large t 
in ( [III ) and keep the multipliers with the number of re- 
fiections n < n* — /3/(2iL7r) only to come to the long 
time asymptotics: 



< gi0c(K,i)g-J0c(s/,O) g7(Tt/T)C_^ 

^tL- 



a -.2g, 



' sinh(Tx{t — ia) / (3)sinh{Tr{—t + ia) / j3) 



x)iL 



(15) 



where zITl/T) = r^/(*^^) and i{Tl/T) approach the 
constant "f{oo) on the order of 1 as \n{TL/T)z{TL/T). 
Our asymptotic analysis following in essential Maslov's 
paper |10|| shows that the low energy exponents approach 
their free electron values as exp[Ti^ lnr/(T7r)]. The effect 
accounts for prolongation of the paths due to the finite 
reflection. In particular, it determines the coefficient c{g) 



of the corrections to the non-universal zero tempera- 
ture value of the conductance variation due to impurities: 
AG.mp (X -iL/l)iT/EFy-H'^ - c{g){pTL)'^) in a uni- 
versal way Q . After substitution of ( p^ into Eq. (|l^) , 
the current suppression produced by the even to terms 
of the interaction becomes equal to: 



< A^m > = 



r(TO2(i-z)) 1^7 



i?,„2, {2mkmFL)TL 



ttT 

t7 



Ef 



'(9-1) 



(16) 



where function fa{x) ~ sinh(a;) J|_|_ r(a ± ix/n) characterizes the V — T cross over. It approaches T^{a)x{l 
(lnr(a))'(x/7r)2) at a; ^ 1 and 7r(x/7r)2°~^ at x 1. Function R specifies the kF — \/L crossover as: 



R2g{x) 



nV\l + 2rg) . 



'^i/2+2r<,(a^/2) - r2(l + 2rg) | j//^t(^/2 - ^ > 1 ^^^^ 



It brings out an interference structure in the conductance 
versus the chemical potential at low energy. This struc- 
ture coincides with the one of Eq. (^ at to = 1. At 
larger to, however, the oscillations are more frequent. In 
particular, there can be the unchanged L/a number of 
maximums of the conductance in between its neighbor 
minimums at the half-filling and the integer filling. This 



effect has not been accounted for in . 
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FIG. 2. Schematic temperature dependence of the conduc- 
tance corrections produced by the m — 2 Umklapp interaction 
AG2 at (7 < 3/4 (linesl) and at g > 3/4 (hnes2) and by the 
random impurity potential AGimp [pII1C|] (lines 3) . Solid lines, 
Ethr ~ 0; dashed lines, Ethr ^ Tl- The dot lines 1 and 2 (1' 
and 2') are the full conductance correction at Ethr ~ (finite 

Ethr)- 



The odd m terms of the spin electron current will 
meet Eq.(p^ after substitution + 1 instead of rn^ 
into the powers and the index of the / function in 
this equation. In the spinlcss case we have to sub- 
stitute 2m^ there and into the index of R. Com- 
bining the above results we can outline a tempera- 
ture dependence of the conductance correction produced 
by the m = 2 Umklapp interaction (Fig.^). For 
spin electrons its magnitude increases/decreases follow- 
ing {Ef/Tl){T/Ef)'^^'^ as T going down above Tl 
and follows {T /TLf{TL/ EfY<)-^ , if Ettr < Tl; oth- 
erwise, the correction starts to decrease exponentially 
exp{—Ethr/T) below Ethr and keeps on decreasing like 
{T/TL)HTL/EFy^-HTL/Ethr)^+'''-<>sm^{2k2FL - nrg) 
below Tl- The T(> Tl) dependence is similar to that 
of the conductivity of infinite wire found by Giamarchi 
Similar dependence with T replaced by V could be 
predicted for the zero temperature differential conduc- 
tance dG2{V) atV <Tl- 

In summary under perturbative condition we have de- 
scribed a hierarchy of the threshold features produced 
by the Umklapp backscatterings at the rational values of 
the occupation number inside the ID channel connecting 
two Fermi liquid reservoirs. In the differential conduc- 
tance (its derivative) vs. the chemical potential /thresh- 
old energy at a finite voltage, the threshold structure is 
an asymmetric peak of width max{T,TL) located at the 
crossover Ethr ~ as the chemical potential moves away 
from the rational filling. These peaks are produced by 
any repulsive interaction at the half-filling of spin elec- 
trons in the wire. In the conductance vs. temperature, 
we predicted a maximum below Ethr due to crossover 
from the Umklapp backscattering to the impurity sup- 
pression and an asymmetric minimum at Ethr if the in- 
teraction is strong enough. However, if the interaction is 



weak so that g > 3/4 the suppression of the conductance 
even at the half-filling may be difficult for observation 
while tii/Tl <C 1. 



IV. ID MOTT-HUBBARD INSULATOR: 
NON-PERTURBATIVE RESULTS. 

In this section, to clear up the difference between the 
Mott-Hubbard insulator and the band gap one, we map 
the problem at low energies and at high temperature 
onto the exactly solvable models making use of a free 
fermion value of the constant g of the forward scatter- 
ing inside the wire. The results are shown in Figs. 3 
and 4. At low energies when T,V <^ Tl {T: temper- 
ature; V: voltage), we have found that a new energy 
scale Tx oc TLexp[—2^Jrin? — /Tl] appears in the sys- 
tem if the chemical potential /i = Ethr/i'^g) of the wire 
is small enough: m? — IJ? /Tl ^ 1. Below Tx the con- 
ductance is not suppressed and the current increases lin- 
early. Above this energy the current saturates and the 
conductance goes down as T^/T reaching small values 
K. exp[-2y/rn?^^ /Tl] aX. T k Tl- At high tempera- 
ture T ^ TL,m we confirmed the asymptotical behavior 
of the conductance: G= (1 — cst^'ll+cosh ^)~-^)/7r for a 
Mott-Hubbard insulator which has been found in the pre- 
vious section in the perturbative regime of a small gap 
|l8t . A brief physical explanation to these results fol- 
lows. At low energies T < Tl and /i <C m the charge 
field is quantized inside the wire at its values related 
to the degenerate sin-Gordon vacua. Rare low energy 
excitations tunnel through the wire with the amplitude 
cx exp[—m/TL] as (anti)solitons switching the quantized 
value of the field. The whole process of tunneling, how- 
ever, includes transformation of the reservoir electron 
into the sin-Gordon quasiparticles and back. This trans- 
formation results in a non-trivial scaling dimension of the 
tunneling operator equal to 1/2 for the Mott-Hubbard in- 
sulator connected to the Fermi liquid reservoirs indepen- 
dently of any parameters. In the case of the band insula- 
tor, this dimension is marginal (= 1): the transformation 
is trivial and does not introduce additional energy depen- 
dence. The infrared relevantness of the tunneling with 
the 1/2 dimension brings out above resonance at zero en- 
ergy. Meanwhile, the exponentially small tunneling am- 
plitude specifies the narrow width of this resonance equal 
to the crossover energy. Increase of favors tunneling 
of the quasiparticles of the same sort and ultimately pro- 
duces their finite density in the wire. Then the inter- 
action between these quasiparticles described with the 
two-particle 5'-matrices |Q dependent on g emerges. At 
low momenta the 5'-matrix for the quasiparticles of the 
same sort is inevitably free fermion like, as at g = 1/2. It 
manifests in the renormalization group (RG) flow derived 
from the Bethe anzats solution for the massive phase p^ ] 
of the sin-Gordon model and in the exponent calculated 
for the Tomonaga Luttinger liquid (TLL) phase at low 
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density ||2^. Increase of T, on the other hand, is expected 
to entail, first, a thermaUy activated behavior of the con- 
ductance (X exp[—2m/T] a.t Tl < T < m ^ and then 
a power law dependence at m < T. Since the effective 
value of g, in general, scales with energy, the 1/T depen- 
dence we found for g = 1/2 may vary at higher energies 
T ^ m depending on the high energy value of g. 

Transport through the finite length wire under a 
constant voltage V between the left and right leads 
could be described in the inhomogeneous Tomonaga- 
Luttinger liquid model (TLL) with the Lagrangian 
[||,|3|: / dx{Y,b^b{x,(pb,dt(l)b)+ Cts{x,Vt,(j)c,(j)s)} as- 
sociated to the Hamiltonian (0,0). The bosonic fields 



J dxCt 

where Lp{x) ~ 6{x)6{L — x) specifies a one channel 
wire of the length L adiabatically attached to the leads 
X > L,x < Q and vf{Ef) denotes the Fermi veloc- 
ity(energy) in the channel. The parameter Ak2F varies 
the chemical potential ^ of the wire from its zero value 
at half-filling. In a real experiment as we discussed before 
this chemical potential is linearly changed by variation of 
the electrochemical potential of the screening gate or the 
average of the reservoir potentials. Outside the Hubbard 
gap this parameter coincides with the momentum trans- 
ferred by the backscattering: four Fermi momenta minus 
a vector of the reciprocal lattice, and relates the present 
results to the ones of the previous section. We as- 
sume > below. The constant of the forward scatter- 
ing varies from gdx) = g inside the wire {x € [0, L]) to 
gc{x) = Qoo = 1 inside the leads, and the Umklapp scat- 
tering of the strength U is introduced inside the wire. 
The charge velocity Vc{x) changes from vp outside the 
wire to a some constant Vc inside it. In the absence of 
the Umklapp scattering, Vc — vpjg and < g < 1 is de- 
termined by the forward scattering amplitude of the bare 
short range interaction between electrons. Approaching 
the half-filling put the Umklapp scattering on. It entails 
an essential renormalization of the low energy value of 
which flows to its free fcrmion value g = 1/2 in the 
massive phase (/i < to) where the coefficient of the 
cos-term scales to ~ rv? and on approaching this phase 
p6[ \ TO. This value of g will be assumed below. 
The zero frequency current through the wire equals / — 
V/t:+ < ibsc >, where the backscattering current ||2^ is 
ibsc = -2ElU/vF Jq dxsm{AkFX + 2Vt + V2(j)c{x)). It 
will be shown later that 2'kEfU is a half gap to, opened 
by the backscattering (Q) in the charge mode spectrum 
inside the wire. 

l.High temperatures T > TF,m - The average backscat- 
tering current < hsc >= J D(f>Ibscexp{i J dt J dx{Cc + 



(j)b{x,t),b — c, s relate to the deviations of the charge 
and spin densities from their average values as follow- 
ing: pb{x,t) = {dx4'b{x,t))/{V2TT), respectively. The 
first part of the Lagrangian describes a free electron mo- 
tion modified by the forward scattering interaction. The 
second part of the Lagrangian introduces backscatter- 
ing inside the wire. Only its term corresponding to the 
Umklapp process of four Fermi momenta transfer is im- 
portant near half-filling. This term does not involve the 
spin field. Therefore, our consideration will be restricted 
to the charge field only. For the clean wire this field is 
characterized by the Lagrangian: 



= y ^-[^{^(^l^j - [^t^) }-^^i-)cos{Ak,FX + 2Vt + V2Mt,m (18) 



^bsc)} can be written as a formal infinite se- 
ries in U. Each term of it is an integral 
of product of the free bosonic correlators < 
exp{i^/2(f){x,t)}exp{—iV2(l}{y,0)} >. Such a corre- 
lator approaches its uniform TLL expression when 
X, L — x,y, L — y ^ Vc/T. Substitution of this form into 
the above series allowed us Q to find L-proportional 
part of the backscattering current neglecting the bound- 
ary contribution in the perturbative case. However, 
the problem is not perturbative, in general, due to a 
finite gap 2to creation. Therefore, application of the 
uniform correlator will give us a part of the backscat- 
tering current cx min(L, Vc/m) with the relative error 
0{max{TF,m)/T), which is of the order of ratio of the 
border piece oc Vc/T to the essential part of the "bulk" 
one. This relates to the high-temperature asymptotics 
of the whole current. 



7Z G 




Xx Xl Km/2 X 

FIG. 3. Schematic linear bias conductance G vs. tempera- 
ture near half-filling fi <^Tl: curve 1 for the weak interaction 
m ^Tl, curve 2 for the strong one. 



Calculation of above series with the uniform TLL cor- 
relator is equivalent to expanding the value g — 1/2 into 
the leads. Following Luther and Emery we map this 
bosonic Lagrangian onto the free massive fermion one 
p9|,p0[ with the density of Lagrangian 
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Hi 



a=R,L=± 



(19) 



right (left) chiral fermion 
backscattering current I^sc 



field. 



Here i/'-R(l) is 
The fermionized 
2mi dx[exp{i{4:kpx + 2Vt)}il}'l^{x,t)'ipji{x,t) — h.c] 
is the doubled backscattering current for the fermions 
p3| under doubled voltage. To find its average we just 
need to know the fermionic reflection coefficient i? as a 
function of dimensionless energy uj = e/m: 



R{lu) 



sin^(Vw2 _ it^) 



{uj^ - 1) + sin^yAJ^^tL) 



(20) 



where = m/Tj^ denotes the dimensionless traver- 
sal time. The analytical continuation is assumed for 
< 1. Since the chemical potential for the right/left 
chiral fermions is fj, ± V, respectively, the total current 
can be expressed as 

V m sinh(WT) 
I = — — - X 

TT TT 



duj 



R{uj) 



cosh((mw - Ai)/T) + cosh(V"/T) 



(21) 



Only the leading term in max{m, Ti,)/T of the right hand 
side of (^) is meaningful. Extracting it, we find the high- 
temperature asymptotics as following 



smh{V/T)B{rntL) 



^ V m 
^ 7 ^ Vcosh((/i)/T) + cosh{V/T) 



(22) 



B{x) 



dhJ- 



sm' 



'■{yUj^^x) 



(^2 - 1) -h sin^(v'a;2 - la;) 



where function B{x) increases as ttx at small x > 
and approaches the constant ~ n at x 3> 
1. Accuracy of this calculation of ( [2l| ) may be 
written as a factor 1 + 0{max{m,T]^)/'max{T,V)) 
to (|ll]) if ImI < max{T,V) or as 1 + 
0{[max{m, Tl ) /max{T, V)] {max{T, V) /Ai)2el''l/^ , oth- 
erwise. The high-temperature conductance (Fig.0) 



TT V T 1 +cosh(^/T) J 



(23) 



approaches zero at T « to if the gap is large enough 
m/TL ^ 1 and < m. This asymptotics gives an ad- 
equate description of the whole high temperature region 
T > m while a high energy value of g was about 1/2. 
Then it complies with the high temperature conductance 
( |l3| ) we have found in the previous section. At others 
values of g we would expect the conductance changing 
its behavior from (|l^) to (^) as g is scaling with low- 
ering the temperature. This crossover could result in a 
non- monotonous dependence of the conductance on the 
temperature, since at high energy g > 3/4 we have seen 
that the conductance is increasing with decrease of the 
temperature. 

2. Low energies T,V << - To find a low energy model 
for our problem we have to integrate out all high energy 
modes. We will try to escape direct integration follow- 
ing Wiegmann's effective way of constructing the Bethe- 
ansatz solvable models [|9| for the Kondo problem and 
for the screening of a resonant level ||2^,Q. First, let us 
substitute 0/\/2 instead of (j> in (^8|). It makes fermions 
interacting inside the leads and non-interacting inside the 
wire. Their passage through the wire may be described 
with the one-electron ^-matrix dependent on the elec- 
tron momentum. The interaction between electrons in 
the leads with some two-particle ^-matrix. Then the so- 
lution could be constructed if the proper commutation 
relations between the 5-matrices are met. Being inter- 
ested in the variation of energy less then around the 
Fermi level, one can simplify the solution keeping the 
one-particle 5'-matrix constant equal to its value on the 
Fermi level. It leads us to the problem of one impurity 
in the TLL. 

For the weak backscattering, the Lagrangian of this 
problem can be written as 



dxCt 



2 %| 



dt(t>c{t-,x) 



An 



dx(l>c{t,x) 



YTlu 



ttvf 



coa{2Vt + V2(t>c{t,0)) 



(24) 



where we rescaled (j) back and introduced a new energy 
cut-off parameter YT^ with dimensionless constant Y 
which will be specified later. Parameter u is related to 
the weak refiection coefficient as:' 



— VpR(^/m). For 
the strong backscattering the tunneling Hamiltonian ap- 



JdxC,=J 



proach may be applied |^l[]. It was associated to the 
dual representation using the field 6 mutually conjugated 
to if) : [0a-{x),(l)a{y)] = i2i:sgn{x — y). The appropriate 
Lagrangian reads 



2 ^vl 



dtOc{t,x) 



Att 



dxOc{t,x) 



Att 



YTlu' 

TTVp 



s{vt + es,o)/V2) 



(25) 
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with u'^ — v\.[\ — R{iJL/m)) proportional to the free mas- 
sive fermion transmittance and the vohage multipHed by 
g factor Both these Lagrangian are, indeed, equiva- 
lent js^] if interaction dependent relation between u and 
u' is met [|4|j3|. The above model (|2|) or@ charac- 
terizes the point scatterer of any backscattering strength 
at low energy |Q. Although, the exact relation between 
u or u' and the bare parameters of the scatterer remains 
unknown. Our problem is dually symmetrical to that 
of Kane and Fisher: suppression of the direct current 
in their problem equals suppression of the backscatter- 
ing one in our case. This correspondence allows us to 
re- write their solution IsljTil as follows: 



T 

I ^ — arctan(F/T:^,), T = 

TT 

T 1 T 



(26) 



where ip denotes the digamma function and satisfies: 
'ip'{l/2) — 7r^/2, Tp'i^) oc 1/2^, X — + cxD, and a new energy 
scale Tx varies from = YTl/ VaR at the weak 
backscattering (|^ to = YTl{1 - R)/t^ at the strong 
one @. 

Let us, first, compare this result with the perturbative 
one H^] of the previous section. The latter was derived 
making use of the long-time asymptotics for the correla- 
tor (|l|): 



< e 



>— const (^aTL)^^ { 



(ttT/Tl 



sinh{'K{t — ia)T)sinh{'K{~t + ia)T) 



)Fix)Fiy) 



(27) 



where a = 1/Ep and F was simplified in the previous 
section as: F{x) = (.t/L)^'- n±,™=i("i ± x/L)^'^""*' « 
const' X [x{L — a;)/i^]^'' and const' — e'^'^°°\ One can 
see that substitution of this asymptotics in the whole 
formal series for the backscattering current discussed 
above implies transformation of the Lagrangian (^8|) into 
the one of ( p^ with the coefficient for the cos-term: 
e'^r{l + r)mJr+i/2{EthrtL)/{V^{2EthrtLy+^^'^} instead 
of il3iii g^jjj]^ another energy cut-off Tr. This model would 
be equivalent to that we constructed before in the weak 
perturbative regime if we can meet 

sin(^ti) V7rr(l + r) Jr+i/2{t^tL) 

At zero r it exactly specifies 1^ as a constant on the or- 
der of 1. However, if r 7^ (r = 1/3 for g = 1/2 we 
assumed), Y increases oc {^t^Y at large /ii^. Moreover, 
there is a mismatching between the oscillating structures 
of Jr+i/iip-th) and sin(/xii) which cannot be naturally 
accounted for by a smooth variation of the energy cut- 
off, but sooner by a small deviation of the traversal time 
ti^ in the free electron reflection coefficient (|2^) from its 
bare value L/vc as changes. Such behavior results from 
penetration of the interaction inside the wire. It is de- 
scribed by the finite reflection coefficient r in the inho- 
mogeneous TLL model. The phenomenon is more im- 
portant for /i 3> max{TL, m) when the electron propaga- 
tion through the wire is not suppressed. In the opposite 
regime of small |/x|/Ti < 1 no interference structure is 
expected and Y remains constant. Finally, under this 
choice ( p8| ) of Y one can see that ^ T^. Therefore, 
the solution ( |2^ ) coincides with the perturbative result 
(|l|, |l^) that is / - F/tt cx -V^ and G - I/tt cx -T^. 



A TCI 




FIG. 4. Schematic current / vs. voltage V near half-filling 
<C Tz,: curve 1 is zero temperature dependence, curve 2 is 
the high temperature T ^ Tl one, the dashed lines are the 
low voltage / = V/iv and high voltage asymptotics. 

Turning to the case m/TL > 1 we cannot use the per- 
turbative expression ( p8| ) anymore: The perturbative se- 
ries is not convergent due to a finite gap creation. Then 
above non-perturbative consideration is necessary. Ap- 
plication of the solution ( p6| ) in this case reveals a quite 
remarkable property of low energy transport through the 
Mott-Hubbard insulator. There is an exponentially small 
value of = (1 - R{fi/m))YTL/7: oc TLexp{~2m/TL) 
for <C ?Ti. Hence, the zero temperature current / (Fig.^) 
is not suppressed for the voltage less than and sat- 
urates at Tx/2 value when < V < T^. Similarly, 
the conductance (Fig.^ displays a small decrease oc 
below its zero temperature value I/tt with increase of 
T for T < Tx and approaches its exponentially small 
asymptotics G = T^/Ut) oc exp{-2m/TL)TL/T above 
Tx < T < Tl. As 1^1 increases, the reflection coefficient 
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R{lJ,/m) on the Fermi level goes down and exceeds 
Ti, finally, approaching its weak backscattering value 
2Tx — YTl/ y^R{fi/m), where the perturbative consider- 
ation is applicable. 

In summary, we studied transport through a ID Mott- 
Hubbard insulator beyond perturbative approach. As- 
suming that g = 1/2 near the half-filling in agreement 
with the Bethe ansatz solutions we mapped the prob- 
lem onto the exactly solvable models and found cur- 
rent vs. voltage at high temperature T > max(m,TL) 
and at low energy T,V < Tl- The solution of these 
models shows, in particularly, that the high-temperature 
transport through the Mott-Hubbard insulator is simi- 
lar to the one through the band gap insulator aX, g — 
1/2. At low energies, however, there is always a regime 
where the transport remains non-suppressed in the ab- 
sence of the impurity backscattering. For the strong 
interaction resulting in the opening of the large Mott- 
Hubbard gap, the transport through the wire is sup- 
pressed near the half-filling almost everywhere inside the 
gap except for an exponentially small low energy region 
V,T <TLexp{-m/TL). 
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